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Spin fluctuations and superconductivity in layered /-electron superlattices 
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We investigate magnetic and superconducting properties of layered /-electron superlattices within 
the fluctuation exchange approximation (FLEX). We show that spin fluctuations, which are char¬ 
acterized by the maximum value of the spin susceptibility in the 3-dimensional (3D) Brillouin zone, 
are strongly suppressed in /-electron super lattices. However, effective 2D spin fluctuations can 
be increased due to the spatial confinement of the /-electrons. Therefore, the tendency towards 
da, 2 _^ 2 -wave superconductivity, mediated by these spin fluctuations, can be strongly increased in 
/-electron-superlattices. This is in sharp contrast to superlattices composed of conventional s-wave 
superconductors, where superconductivity is generally suppressed. 

PACS numbers: 


I. INTRODUCTION 

Recent experimental realizations of layered superlat¬ 
tices, Celns/Lalns and CeCoIns/YbConIns, have opened 
new possibilities in the field of /-electron systems.^^l® 
Due to a non-trivial interplay of strong correlations and 
tunable dimensionality, novel phenomena have been ob¬ 
served in these /-electron superlattices which have not 
been seen so far in existing magnetic/superconducting 
superlattices composed of weakly or non-interacting sys¬ 
tems. EHm 

For example, it has been found that magnetic proper¬ 
ties of Celn3(n)/Laln3(4) superlatticeJ^l depend on the 
thickness of the Celn 3 -layer within the unit cell of the 
superlattice. In bulk Celn 3 , the coherence temperature 
is Tcoh ^ 50(K) and the Neel temperature is Tat — 10(K) 
with an ordering vector Q = (tt, tt, Remarkably, 

it has been reported that the Neel temperature of the 
superlattice is suppressed when the width of the Celn 3 - 
layers, n, is reduced, and eventually approaches zero for 
n = 2. At the same time, a linear temperature depen¬ 
dence in the in-plane resistivity is observed for n = 2, 
suggesting that the dimensionality of the antiferromag¬ 
netic (AF) spin fluctuations is reduced from 3-dimensions 
in bulk Celn 3 to 2-dimensions in the superlattice. Such 
an anomalous behavior in the resistivity has never been 
observed in previous studies of magnetic superlattices, 
and would be characteristic for /-electron superlattices. 

Furthermore, superconductivity in CeCoIns/YbCoIns 
superlattices has been investigated.bulk CeCoIns, 
the coherence temperature is Tcoh ^ 50(K) and 

the su perco nducting transition temperature is Tc ^ 
2.3(K).™llCeCoIn5 is located near an AF quantum crit¬ 
ical point and AF spin fluctuations are expected to be 
important for the normal state as well as for supercon¬ 
ductivity. In the bulk system, the AF spin fluctuations 
are especially strong around Q = (7r,7r,7r) due to the 
nesting of the Fermi surface, and they can be character¬ 
ized as 3D-like.lI^J^^ It is generally considered that the 
superconductivity exhibits symmetry and is 

mediated by these AF spin fluctuations. Experiments on 
CeCoIns /YbCoIns superlattices have demonstrated that 


superconductivity exists even for thin CeCoIns-layers 
and that the superconducting transition temperature, Tc, 
is suppressed as the width of the CeCoIns-layers in the 
unit cell is reduced. However, it must be noted that at the 
same time effects of disorder, which are estimated from 
the residual resistivity, are increased in thin CeCoIns- 
layers. It is thus unclear, how the superconductivity be¬ 
haves in “clean” /-electron superlattices. 

Motivated by these experiments, there have been sev¬ 
eral theoretical studies. The effects of a possible Rashba- 
like spin-orbit coupling due to local inversion symme¬ 
try breaking near the interfaces of the Ce-layers and 
the spacer layer^^^H^ have been investigated. When 
the Rashba-like interaction is sufficiently large, the Pauli 
depairing effect is greatly suppressed and novel super¬ 
conducting states might be stabilized when a magnetic 
field is applied. In another theoretical study, the ex¬ 
perimental data was analyzed based on the Berezinskii- 
Kosterlitz-Thouless transition by regarding the superlat¬ 
tice as a junction composed of a normal metal and a 
superconductor.!^ If this junction picture is applicable 
to the /-electron superlattice, then superconductivity in 
the YbCoIns-layer would be strongly suppressed, because 
of a large mismatch between the Fermi velocities of the 
CeCoIns-layer and the YbCoIns-layer, leading to a 2- 
dimensional superconductivity in the CeCoIns-layers. 

In these previous studies neither electron correlations 
nor the superlattice structure are explicitly considered. 
However, these are two key ingredients in /-electron 
superlattices and distinguish them from all the exist¬ 
ing non-interacting superlattices and the bulk /-electron 
compounds. To understand /-electron superlattices, it 
is necessary to clarify the impact of electron correlations 
and the superlattice structure, and also their possible 
interplay. In two previous studies, the present authors 
already d iscus sed the Kondo effect and quasi-particles 
propertie d^^ * ^^ ! using the dynamical mean field theory 
(DMFT), which captures local strong correlations, but 
neglects non-local fluctuations. This time, we analyze 
magnetic and superconducting properties of /-electron 
superlattices using the fluctuation exchange approxima¬ 
tion (FLEX) in order to describe spatially extended spin 
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fluctuations!^ We use a periodic Anderson model (PAM) 
which is deflned on a superlattice. This model can be con¬ 
sidered as a minimal model to describe /-electron super¬ 
lattices, because it takes into account both the electron 
correlations and the superlattice structures. 

This paper is organized as follows: In Sec. [Ilj we in¬ 
troduce our model, the FLEX approximation for the spin 
fluctuations, and the Eliashberg equations for the super¬ 
conductivity. Spin fluctuations are discussed in Sec. Ing 
and the superconducting instability is examined in Sec. 
m Einally, in Sec. |V|we shortly summarize this paper. 


II. MODEL 


Because the /-electrons in the Yb-sites form a closed 
shell, the YbCoIns-layers in CeCoIns/YbCoIns super¬ 
lattices can be treated as normal uncorrelated metals. 
Indeed, the resistivity in bulk YbCoIns shows a mono¬ 
tonic temperature dependence without any signature of 
the Kondo effect. Similarly, in the Celns/Lalns super¬ 
lattices, only the Celns-layers provide /-electrons at the 
Eermi energy. Therefore, both superlattices can be con¬ 
sidered as heterostructure composed of layers including 
/-electrons and layers without. In order to understand 
these /-electron superlattices, we introduce a PAM which 
consists of two kinds of layers, henceforth called “A- 
layers” and “B-layers”. The A-layers include conduc¬ 
tion electrons (c-electrons) and /-electrons, which cor¬ 
responds to Celns- or CeCoIns-layers, while within the 
B-layers only c-electrons exist, which corresponds to the 
Lalns- and YbCoIns-layers. Our Hamiltonian thus reads 


H = 


rr'a 




Y ^ivflafr'c 

rr' EA,a 


ra rEA,a 

T ^ [^rafrcT + fra^ra] + U 

tGAjCt rGA 


( 1 ) 


where Cro- and fra are annihilation operators for the 
conduction electrons and the /-electrons, respectively. 
r = = {x^y^z) is a site index which is composed 

of an in-plane index , ry, and a layer index z. a corre¬ 
sponds to the spin index. Each layer forms a square lat¬ 
tice, and hopping is only allowed between nearest neigh¬ 
bor sites for simplicity; ^ for ry 7 ^ Cy and 

'^r\\zr\\z' ~ ^2 z ^ z' where a = c^f. The number of 
A-layers and B-layers within the unit cell of the super¬ 
lattice are given by La and Lb- In the present study, 
we flx La = 1 for which effects of the spatial conflne- 
ment of the /-electrons are expected to be particularly 
strong. Because of ^2 = 0 and the absence of a direct 
hopping between the A-layers, which are separated by the 
B-layers, the /-electrons can move along the ^-direction 
only through the B-layers. The model parameters are 
chosen as (^i, E, 5 /) = (5.0, 2.0,0) and the total fllling is 


flxed to n = = 0.95, which are a reasonable set 

of parameters and similar to the ones used in the previ¬ 
ous DMET study.l^ The interaction strength is flxed at 
a moderate value, U = 3.0, for which a clear divergence 
in the spin susceptibility of the 3D system {Lb =0) is 
visible at low temperature. The 2 :-axis hopping char¬ 
acterizes an anisotropy of the system for = 0 . 

In order to analyze momentum-resolved properties, we 
perform a Eourier transform. 


^jziZ2,k\\kzl^k\\kzlcr^ 

k\\kzl 

(2) 

^ ^ ^jziZ2,k\\kzf^\\^z^^ 

k\\kz 

(3) 


where the unitary matrices and U’^ are deflned as. 


Tjc 

jziZ 2 ,kukzl 


jziZ 2 ,kukz 


gi/c||r|| ^ikzZ-\-iqfz2 


gi/c||r|| ^ikzZ 


(4) 

(5) 


k\\ = {kx^ky). The layer index 2 : is parametrized as z = 
Lzi + Z 2 with t) < Z 2 < L ioT and Z 2 = I for 
and 0 < / < 1/ for qf = 21:1 jL. The momentum 
along the z-axis is deflned within the reduced Brillouin 
zone (RBZ), 0 < kz < 27:jL. Yy is the total number of 
sites within a layer and Nz is the total number of layers. 
Thus, the total number of sites is given by Y = Yy Y;^. 

In ELEX, we focus on spatially extended spin fluctua¬ 
tions. I^The selfenergy and susceptibilities in the normal 
state are given by. 


j:^{k) 

Vf{q) 

Xoil) 

xi{q) 

xiiq) 


^Y.Vf{q)Gf{k-q), 

q 

(6) 

\u^xl{q) + ^u^xiiq) - u^xliq) 

(7) 

q 

(8) 

Xo(9) 

l + UxliqY 

(9) 

Xo(9) 

i-UxkqY 

(10) 


where k = {iuj^^k) with k = {k\\,kz). {k) is the /- 

electron Green’s function in the //c^-basis. Note that the 
/-electron contributions to the total spin susceptibility 
are dominant, especially near magnetic criticality, and 
that they are strongly enhanced by the interaction U. 
On the other hand, contributions from the c-electrons 
are not enhanced by an interaction term in the present 
model. 

The superconducting instability is investigated within 
the linearized Eliashberg equation for the singlet gap 
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FIG. 1: Spin susceptibility 


0, q) at T = 0.1 in the xz-plane at ^2 = ti foi* different Lb- 


function (k)^ 

= -^ E - k')\Gk{k')\^Ak{k'), (11) 

k' 

V/iq) = U- lu\liq) + lu\f{q). (12) 

It is noted that, c-electrons can only become supercon¬ 
ducting via the /-electrons by the proximity effect, be¬ 
cause of the absence of c-electrons interactions. The 
proximity effect is well taken into account in our cal¬ 
culations because the /-electron Green’s function G-^ 
fully includes the hybridization processes between the /- 
electrons and the c-electrons through V. 


III. SPIN FLUCTUATIONS 

In this section, we discuss the spin fluctuations as cal¬ 
culated by the FLEX. First, we consider an isotropic pa¬ 
rameter set, ^2 = tf, where anisotropy between the xy- 
and 2 ;- directions can only originate from the superlattice 
structure when > 1. Figurel^shows the g^-dependence 
of the magnetic susceptibility, at rela¬ 

tively high temperature, T = 0.1. The positions of the 
maximum values of oscillate depending on the num¬ 
ber of spacer layers Lb- When is even, the maxi¬ 
mum values are located at Q = (tt, tt, tt/L); when is 
odd they are at Q = (7r,7r,0). These momenta corre¬ 
spond to the spin configurations shown in Fig. When 
is odd (even), the magnetic coupling between dif¬ 
ferent A-layers, which is mediated by the c-electrons, is 




FIG. 2: Schematic picture of the spin configurations for 

Ls = 1 (left panel) and Lb —2 (right panel). Shaded layers 
and white layers are A-layers and B-layers, respectively. 


ferromagnetic (antiferromagnetic). Due to the proxim¬ 
ity to the A-layers, small moments are induced into the 
spacer B-layers in a consistent way with the magnetic 
structures of the A-layers. Similar oscillating inter-layer 
magnetic structures and induced moments have been 
found in DMFT calculations, which support the present 
FLEX study.l^Eurthermore, such oscillations in the mag¬ 
netic inter-layer coupling have also been commonly found 
in ferromagnetic superlatticesOne intuitive under¬ 
standing of this phenomena is based on the RKKY in¬ 
teraction bet ween magnetic layers separated by metallic 
spacer layers.^^^^^^ The magnetic inter-layer coupling is 
asymptotically given by ^ J^\ii2kFZ /with the Eermi 
wavenumber kp and coupling strength J. In the present 
study, the system is close to half filling so that the Eermi 
wavenumber of the c-electrons at U = 0 along the 2 )-axis 
is ^ 7r/2, which leads to the above-mentioned periodicity 
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FIG. 3: Temperature dependence of the maximum value 
of = 0,q) for t 2 = t\. Shaded region indicates the 

characteristic temperature scale Tq. 


in xiio)- 

We can estimate the strength of the spin fluctuations 
by the maximum value which we show 

in Fig. We want to note here, that the magnitude 
of this maximum value is strongly parameter dependent, 
because already small changes in the density of states 
can lead to a substantial enhancement of the Stoner fac¬ 
tor (Eq. , if the system is close to magnetic criticality. 
A general behavior observed in our calculations is that 
the susceptibility for the 3D bulk system (without su¬ 
perlattice structure, = 0 ) strongly increases below a 
characteristic temperature Tq, indicated by the shaded 
region in the figure, Tq ^ 0.05. At the same temper¬ 
ature, the single peak in which is present at high 
temperature (Fig. [^, is split as shown in Fig. Due to 
the hybridization between c- and /-electrons, the Fermi 
surface is split and the nesting properties are changed 
at low temperatures as seen in Fig. Thus, the spin 
fluctuations are strongly affected by E for T < Tq where 
heavy fermions with long lifetime are well formed within 
the present FLEX calculations. 

In the /-electron superlattice, 7 ^ 0, the Fermi sur¬ 
face differs strongly from the 3D bulk system.!^ Par¬ 
ticularly, a g'; 2 -dependence of xlio) arises only from 
through V, because a direct hopping between different 
A-layers separated by the spacer B-layers is forbidden in 
the present model. A general behavior observed for the 
superlattice is that, similar to the bulk system, the sin¬ 
gle peak in the susceptibility at high temperature is split 
into four peaks at low temperatures, T < Tq ^ 0.05, as 
shown in Fig. However, the abrupt increase in the 
susceptibility, which is observed in the bulk system at 
low temperature, is cut off in the superlattice around Tq 
and the susceptibility decreases for T < Tq, as shown in 
Fig-i While at high temperature, T > Tq, the Fermi 
surface is mainly determined by the c-electrons, below 
To the /-electrons and thus the superlattice structure 
become important. As a consequence, for temperatures 
below To the nesting of the Fermi surface in the super¬ 
lattice is changed and Q ^ (tt, tt, tt) is no longer a good 
nesting vector as can be seen in Fig.[^ Therefore, as T is 


decreased and the Fermi surface becomes affected by the 
superlattice structure through V, the lack of good nesting 
properties of the Fermi surface cuts off the enhancement 
of the spin fluctuations. We note that the large value of 
max[x{] for T^ = 1 originates in a large density of states 
near the Fermi energy due to the superlattice structur^^ 
which leads to a substantial enhancement in the Stoner 
factor, as mentioned above. However, while this strong 
enhancement at finite temperature for T^ = 1 is strongly 
parameter dependent, the decrease of the susceptibility 
in the superlattice below Tq has been observed for a wide 
range of the parameters. 

As shown in Fig. spin fluctuations in the superlat¬ 
tice are smeared out at low temperatures within the Qz 
direction. The -dependence of x{(^) becomes weaker 
as Tb is increased. Thus, spin fluctuations become more 
2-dimensional-like when T^ is increased. However, this 
also means that the peak height of (g) in the 3D Bril- 
louin zone is not an appropriate measure for the strength 
of the spin fluctuations when T^ is large. In order to es¬ 
timate the strength of the spin fluctuations, we consider 
an effective 2D spin susceptibility 

xi2D(.^<^n,q\\) = (13) 

x{ 2 d(^||) maximum at Qy ^ ( 7 r, 7 r) for any Lb in 

the present model. We show the temperature dependence 
of the maximum values of xf 2 D(^l|) Although 

max[x{ 2 D] superlattice (T^ > 1 ) does not show a 

divergence at low temperature, this effective 2D spin sus¬ 
ceptibility is clearly enhanced compared to the 3D PAM 
{Lb = 0) for the calculated temperature range. We note 
that the change in dimensionality of the spin fluctuations 
from 3D to 2D has also been experimentally observed in 
the Celus/Lalus superlattice as the Celus-layer thick¬ 
ness {La in our model) was tuned with a fixed Lalus- 
layer thickness (T^).® Although we cannot directly com¬ 
pare our results to the experiments, the suppression of 
the magnetic order in the superlattice and the calculated 
2D-like character of the spin fluctuations are consistent 
with the experiments. 

Up to now, we have analyzed the magnetic susceptibil¬ 
ity for an isotropic model. However, CeCoIus exhibits a 
cylindrical Fermi surfac^^. In order to investigate effects 
of an anisotropy in the original 3D model, we consider a 
system with a more 2D-like set of hopping parameters, 
^2 = 0.5tf. Similarly to the isotropic hopping parameter 
set, we observe even-odd oscillations of the peak posi¬ 
tions when T^ is changed, as shown in Fig. We show 
the temperature dependence of the maximum values of 
xi (T) within the 3D Brillouin zone in Fig.|q The suscep¬ 
tibility of the 3D PAM {Lb = 0), max[x^(T)], behaves 
again monotonically and rapidly grows at low tempera¬ 
ture. Furthermore, any increase in max[x{(T)] for the 
superlattice, T^ > 0 , is again cut off at low temper¬ 
atures. However, because of the anisotropic parameter 
set, the bulk system, T^ = 0, already includes strong 2D 
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FIG. 4: Spin susceptibility = 0, <gr) at T = 0.03 in the a:z-plane at = tt for —ti. The inset in each panel shows a 

magnification of the BZ as marked by the yellow square in the main plot. 
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FIG. 5: Fermi surface as given by the spectral function 
(—l/7r)ImG-^(z7rT, fc) at a low temperature T = 0.03 with 
/cy = 0 for Ls = 0 (left panel) and Lb = I (right panel). 


spin fluctuations, and the effective 2D spin susceptibility 
^{20 shown in the paper) is less enhanced in the su¬ 
perlattices. Comparing the results for the two parameter 
sets, ^2 = ti and = 0.5tf, we find that the impact of the 
/-electron confinement to the A-layers is stronger when 
the hopping parameters are more 3-dimensional. This 
is reasonable, because in the limit of decoupled layers, 
^2 ^ 0, the superlattice structure does not play a role 
at all. We have confirmed this tendency by performing 
similar calculations for different hopping parameters. 


IV. SUPERCONDUCTIVITY 



FIG. 6: 
Fig-i 


Temperature dependence of the maximum value of 
= 0,qr||) for ^2 = ^ 1 - Symbols are the same as in 


Finally, we investigate the impact of the increased 2D 
spin fluctuations in the superlattice on the superconduc¬ 
tivity by solving the Eliashberg equation 0- First, in 
order to understand the proximity effect in our model, we 
consider s-wave superconductivity by replacing V/ (q) in 
Eq. ([^ by 

= -UoeiujD - (14) 

where we fix = 1.0 and Uq is tuned so that A = 1 when 
Lb = 0 and T = 0.02. We then solve the Eliashberg 
equation without selfenergy. Thus, electron correlations 
are not included in this calculation. In Eig. we show 
the maximum eigenvalues, max[A], of the Eliashberg 
equation which corresponds to the strength of the super¬ 
conducting instability. The largest eigenvalue. A, rapidly 
decreases as is increased, showing slight oscillations 
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FIG. 7: Spin susceptibility Xsi'l^n — 0,q) at T = 0.1 in the xz-plane at = tt for t 2 = 0.5/^. 



FIG. 8: Temperature dependence of the maximum value of 
= 0,q) for = 0.5/?. 


due to changes in the density of states at the Fermi energy 
which depend on the details of the model parameters. 
Thus, if conventional 5-wave superconductivity, mediated 
by phonons, were realized in CeCoIns/YbCoIns superlat¬ 
tices, one can expect that Tc is strongly decreased in the 
superlattice. Such a suppression of the superconductivity 
has been commonly observed in superlattices composed 
of 5-wave superconductors and normal met als.^^MSl Exper¬ 
imentally, Tc is lower in all the previous conventional su¬ 
perlattices compared to the corresponding bulk systems. 
It has theoretically been shown that the superconducting 
transition temperature Tc is suppressed in an exponential 
way in layered systems when the thickness of the normal 
metal layer Ltv is increased, 



FIG. 9: Maximum eigeuvalue A for the s-wave supercouduc- 
tivity at fixed T = 0.02. 


Tc{Ln) — Tc{0) — JTc tanh(I/7v/^o), (15) 

where STc ^ (Tc(I/Ar)/Tc(oc))(Tc(0) - Tc(oc)) and is 
effective coherence length. This suppression of Tc comes 
from the fact that the pairing interaction which exists 
only in the superconductor layer mediates superconduc¬ 
tivity not only in the superconductor layer, but also in 
the normal metal layer. Since this is a general property of 
the proximity effect, one could naively expect a suppres¬ 
sion of superconductivity also for d-wave states, which 
are mediated by the spin-fluctuations. 

In order to examine this further, we solve the Eliash- 
berg equation with the pairing interaction V/(q) cor¬ 
responding to spin fluctuations and the normal selfen¬ 
ergy for ^2 = tf. Contrary to the 5-wave supercon¬ 
ductivity, as seen in Fig. the maximum eigenvalues 
for the dx 2 _y 2 -w 8 ive superconductivity mediated by the 
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FIG. 10: The maximum eigenvalue of the Eliashberg equa- FIG. 12: The maximum eigenvalue of the Eliashberg equa¬ 
tion for the d^ 2 _y 2 -wSive superconductivity when t 2 — t\. tion for the d^ 2 _y 2 -wSive superconductivity when t 2 = 0.5tJ. 
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FIG. 11: Gap functions A{kx,ky,kz = 0) and A{kx,ky = 
0, kz) at LJn = ttT for Lb = 0,1 when t2 = ti and T = 0.02 
in arbitrary units. 


spin fluctuations are enhanced in the /-electron super¬ 
lattices. Although the eigenvalues, niax[A], do not reach 
unity for the calculated temperature range, these results 
suggest that Tc for dx 2 _y 2 -M^me superconductivity can 
be higher in the superlattice than in the bulk system. 
This strong d-waye superconducting instability in the /- 
electron superlattice can be understood by focusing on 
the effective 2D spin fluctuations discussed in the pre¬ 
vious section. In order to stabilize the dx 2 _y 2 -wdive su¬ 
perconductivity with A{k) oc [coskx — cosky), the Qz- 
dependence in V/(q) is irrelevant and we only need to 


consider the dependence. As exemplified in Fig. 11 


typical profiles of the dx 2 _y 2 -w 8 ive gap functions are in¬ 
deed A{k) ^ {coskx — cosky) and their /c^-dependence 
is weak for any Lb- If we neglect the /c^-dependence in 


|G'f(/c)p, the Eliashberg equation is reduced to 

II u/ 

/C|| 

(16) 

= ( 17 ) 


From this equation, it is clear that the most important 
part of the pairing interaction is determined by x{ 2 d(^||) 
which is enhanced in the superlattice (Fig. 1^. This en¬ 
hancement of the effective pairing interaction can lead to 
an increased Tc, which is a consequence of the interplay 
between strong interaction among the /-electrons and the 
confinement of them within the superlattice structure. 
This is characteristic for the /-electron superlattice. We 
note that a similar enhancement of superconductivity has 
been theoretically found in 3D bulk mod els wh en tuning 
the hopping parameter in the z-direct ior l^^F^ I However, 
in these studies the proximity effect does not play a role. 
In our present study, superconductivity is enhanced as a 
result of a subtle interplay between the proximity effect 
and an increase of spin fluctuations in the superlattice. 

Results for max [A] when = 0.5tf are shown in 
Fig. 12 For this parameter set the enhancement of x{ 2 d 
in the superlattice is weak. Therefore, also the maximum 
eigenvalue of the Eliashberg equation are only slightly in¬ 
creased in the superlattice. Tc in the superlattice would 
be similar to the bulk value when the hopping param¬ 
eters are 2D-like. However, the enhancement of x{ 2 D 
in the superlattice due to the /-electron confinement is 
still important for these parameters. It almost cancels 
the suppression of d-wave superconductivity due to the 
proximity effect. Therefore, Tc in the superlattice can 
remain as high as in the bulk {Lb = 0) even for large 
Lb- 


Experimentally it has been observed that Tc is lower in 
the CeCoIus /YbCoIns superlattice than in bulk CeCoIns 
with Tc — 2.3(K). We think that this can be explained 
by two reasons: First, the Fermi surface of CeCoIns is 
cylindricaPIl and thus would be better described by the 





















anisotropic parameter set in our calculations. Second, 
disorder effects seem to be strong for thin CeCoIus-layer 
superlatticesIn the experiments, the thickness of the 
YbCoIus-layers has been fixed and the number of the 
CeCoIus-layers has been tuned. We expect that if the 
thickness of the YbCoIus-layer is changed with a fixed 
CeCoIns-layer width, the behavior of Tc will deviates 
from conventional normal-metal/superconductor super¬ 
lattices. 


V. SUMMARY 

We have investigated the /-electron superlattice based 
on FLEX. We found that the nature of the spin fluctu¬ 
ations is modified by the superlattice structure and the 
Q-vectors corresponding to the maximum in the suscepti¬ 
bility depend on the width of the spacer layers, similar to 
conventional magnetic superlattices. While the strength 
of the 3D spin fluctuations, characterized by max[x{(g')] 


in the 3D Brillouin zone, is suppressed in the superlat¬ 
tice because good nesting properties of the Fermi surface 
are lost, effective 2D fluctuations are enhanced because 
of reduced dimensionality. These enhanced spin fluctu¬ 
ations can lead to higher Tc in the case of da, 2 _^ 2 -wave 
superconductivity in the superlattice than in the bulk 
compounds, which is in sharp contrast to all the conven¬ 
tional superlattice superconductors. We hope that these 
results will lead to further experiments analyzing Tc in 
clean /-electron superlattices. 
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